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Abstract 



There is a growing interest in the possibility that dark matter could be formed 
of weakly interacting particles with a mass in the 100 GeV - 2 TeV range, and 
supersymmetric particles are favorite candidates. If they constitute the dark 
halo of our Galaxy, their mutual annihilations produce energetic gamma rays 
that could be detected using existing atmospheric Cerenkov techniques. 
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1 Introduction 



Galaxies seem embedded in massive spheroidal dark halos [EL 0, pL 0] . Extensions of the standard model 
of particle physics often predict the existence of a new stable particle, a natural candidate for halo dark 
matter. The residual abundance of a particle after the Big-Bang depends on its mass and interactions, 
and to account for the dark matter puzzle a weakly interacting massive particle (W.I.M.P.) is needed, for 
instance a neutrino-like particle with a mass around 1 TeV. Such particles will be very difficult to detect, 
even by accelerator experiments (LHC-SSC) of next century. Halo particles can be directly detected 
when they cross a bolometer detector || ^| but the flux is low above 50 GeV. They can also be indirectly 
detected through their annihilation products [Q, [§| ||, |n], P|, |l2| : annihilations after their capture in the 
Earth or in the Sun producing neutrinos, or annihilations in the halo producing energetic cosmic rays 
(gamma rays, electrons, protons, antiprotons. . .). Existing models for our galactic halo predict a strong 
annihilation rate near the center of the Galaxy. We focus here on the annihilations into photons, because 
photons conserve their initial direction. Each annihilation produces either primary photons (leading to 
a gamma line signal) or secondary photons from the decay of primary annihilation products. The fluxes 
are low again for a large WIMP mass, but this can be balanced by the large acceptance of detectors 
using the atmospheric Cerenkov technique like the existing Atmospheric Cerenkov Telescope (ACT) of 
the Whipple Observatory With similar instruments, it should be possible to detect for instance the 
leading WIMP candidate, the neutralino predicted by supersymmetry [fl4f . 



2 ACT's and backgrounds 

The ACT detects incoming cosmic rays through the Cerenkov light emitted by the shower of secondary 
particles created when the cosmic ray (photon, electron, proton, heavier nucleus. . .) collides with a 
nucleus in the upper atmosphere jl3| . The image shape analysis affords an efficient gamma versus hadron 
separation. The Cerenkov light covers a disk of about 120 m radius on the Earth surface (an area of 45000 
m 2 ), and the detector must be somewhere in this area to see the cosmic ray shower. Less energetic cosmic 
rays produce less Cerenkov light, and the minimal detectable energy (threshold) is linked to the size of 
the mirrors which collect this light. The best observation conditions correspond to a cosmic ray within 
45 degrees from the zenith. As we shall see, the photon flux from WIMP annihilations is largest from the 
galactic center, which can be seen at zenith only if the observatory is built in the southern hemisphere, 
at a latitude near -29 degrees. Observation of Cerenkov light is only possible at night when there is no 
moon, and so the galactic center can be observed for t b s — 2 x 10 6 s per year. The effective angular 
acceptance AJ7 of present ACT's is about 10~ 3 sr with an angular resolution of 3 x 10~ 6 sr, so that they 
are well suited to the study of localized sources. The usual procedure consists in the comparison of the 
flux from an expected source with neighbouring sky regions (on-off experiment). 

To be definite, we shall use as a reference ACT the Granite set-up at the Whipple Observatory, which 
uses two mirrors of about 10 m diameter, 120 m apart, leading at present to a minimal detectable energy 
around 200 GeV. Increasing the number and/or the size of mirrors may lower this threshold. The two 
telescopes effectively survey 63000 m 2 if they are used independantly, but their effective detection area S 
is only 18000 m 2 when the two telescopes are used in coincidence to give the best energy resolution and 
an increased hadronic rejection. The energy of a gamma cosmic ray is determined from a stereoscopic 
view with a precision a(E)/E ~ 0.1 . Adding other telescopes, about 120 m apart, would increase the 
effective area and therefore lower the detectable cosmic flux. The product SAnt i, s is a measure of the 
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efficiency of the telescope to detect cosmic rays, and SASltobs — 3.6 x 10 11 cm 2 sr s for our reference 
ACT. 

There is a large background to the WIMP annihilation signal that we are looking for. Cosmic ray 
hadrons must be rejected, and experience at the Whipple Observatory shows that a rejection factor of 
10 3 should be reached for two detectors in coincidence. The flux of misidentified hadrons then follows a 
power law Jl3| : 

- ^ N '' - = 6 x 1CT 9 (—^—) photon - like hadrons (s cm 2 sr TeV)" 1 (1) 
at dS di 2 dE \ 1 TeV / 

When a hadron is misidentified as a photon, its energy is underestimated, and the preceding formula 
takes this effect into account. Electrons cannot be distinguished from photons with atmospheric Cerenkov 
detectors, as their showers have the same development as photons in the atmosphere. The electron flux 



is measured at energies up to 1 TeV 15, [L6| . The spectrum also follows a power law : 

dt^fdE = 6 X 10 " (l!v) electrons (s cm 2 sr TeV)- (2) 

The electron background dominates the background of misidentified hadrons up to 1 TeV. The fluxes of 
electrons or hadrons are isotropic. Finally, there is the diffuse photon background itself, supposed to be 
due to the collisions of cosmic ray protons on interstellar gas. The spectrum from the galactic center was 
measured by the COS-B satellite JlT], [l8|, (L9j up to 5 GeV : 



df dS dO dE 



= 0.35 x 10~ y I I photons (scm 2 srTeV) -1 (3) 



This power law is the same as for the hadron spectrum, and we may assume that Equation ^ extends to 
higher energies. This is not in disagreement with upper bounds on the photon flux above 5 GeV jl^, |l6| . 
This flux is smaller than the electron background. 

A signal can be detected even if it is small compared to the background, provided that the background 
is measured with a good precision. In our case, an on source-off source method will be available, allowing 
for a reliable estimation of this background. The noise is the uncertainty on the background, that is the 
square root of the number N^ackgr of all photon-like events (including electrons and misidentified hadrons) 
detected during the experiment time t b s in a given energy bin AE, that we shall take as 3a(E) ~ 0.3E. 
For a detector of acceptance SAQ, this number is : 

f d7V 7 dN e dN h 1 Arn 

Nbackar = S Ati t obs { ^ 1 ! 1 1 K AE (4) 

back9r obs \ dtdSdQdE dtdSdfldE dtdSdQdE J v ' 

We must now compare these numbers to the expected signal flux from halo WIMP annihilations. 



3 WIMP annihilation from the galactic center 

We assume that our Galaxy is embedded in a dark halo made of WIMP's of mass m x which annihilate 
in pairs. The annihilation rate is < crV >haio n't, where < aV >haio is the (thermal average in the halo 
of the) annihilation cross-section a times the relative WIMP velocity V, n x = p x /m x is the dark matter 
number density. At large distance r from the center, the halo density decreases as : 

2 i 2 

f \ a +r © (K\ 
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where p© = 0.3 ± 0.1 GeV/cm 3 is the dark matter density in the solar neighbourhood, a — 5 ± 3 kpc is 
the core radius and r = 8.5 kpc is the distance of the Sun to the galactic center EG, EH 



23 . The 



annihilation rate of this dark matter far away from the galactic center turns out to be too low to be 
detected with present size ACT's. 

The dark matter density may however be much larger at the center, as there seems to be a massive 
nucleus at the center of the Galaxy |24 , 25 . According to the isothermal model of Ipser and Sikivie |^6| , 
this nucleus induces an enhancement of the dark matter density within a 150 pc radius around the center, 
with a density at center of : 

p ~ 0.34 a TeV/cm 3 (6) 

The a parameter (0.1 < a < 3) reflects the uncertainties on the rotation velocity of the galaxy, and on 
the contribution of luminous matter to this velocity. Annihilation photons in the TeV energy range are 
unlikely to be absorbed by baryonic matter in the galactic center since, in the worst case of a completely 
gaseous core, a 7-ray crosses less than one radiation length. As seen from Earth, the central region covers 
Ail = 10~ 3 sr, nicely corresponding to the acceptance of present ACT's. We can perform an on-off 
experiment, aiming at the galactic center and comparing with the result off center. The 10~ 6 sr ACT 
angular resolution is small enough to distinguish point sources inside these 10 -3 sr from the more diffuse 
WIMP annihilation cloud. 

Photons appear in the decay chain of the primary annihilation products (quarks, leptons, Higgs scalars, 
W or Z°'s, etc.), and their energy is then continuously distributed. Their distribution dn 1 {E)/dE can 
be computed through a Lund Monte-Carlo (2^] (Figure 1), and depends only weakly on the WIMP model. 
The integral J dn 7 (E) is the photon multiplicity per annihilation. The number of photons due to WIMP 
annihilations within 150 pc from the center, and received on our reference telescope is : 



N^E) dE = St obs dn 7 {E) 



< crV >halo 
47rr 



150 pc 



p(r) 



Airr dr 



(7) 



WIMP annihilations outside the 150 pc sphere, but inside the 10 -3 sr cone, do not contribute to a signal 
defined by an on-source minus off-source subtraction. The cross-sections a and the corresponding energy 
spectra dn 1 (E) / dE , reduced to a delta function in the case of monochromatic photons, are summed 
over all relevant annihilation channels. 

Two different signals can be looked for at the same time. We can look for a gamma line due to the 
direct annihilation of WIMP's into photons, or we can look for the total excess of soft photons due to 
the decays of all annihilation products. 

Soft photons 

The Lund Monte-Carlo codes |27| typically indicate that 50 photons or so are produced in the decay 
chains, but that most of them are at a very low energy compared to the WIMP mass m x (see Figure 1). 
The signal is the number of photons received above the ACT energy threshold E t hresh ■ Table 1 shows 
the integral J dn~/(E) from Ethresh to m x as a function of m x for a quark-antiquark channel (the result 
being very similar for other channels). 



m x (TeV) 


0.4 


0.8 


1.0 


2.0 


4.0 


Ethresh = 0.1 TeV 


0.02 


0.17 


0.32 


1.45 


4.40 


Ethresh = 0.2, TeV 


0.0004 


0.02 


0.05 


0.35 


1.55 
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Table 1 : Number of soft decay photons above threshold. The first line gives the mass m x of the 
particles which annihilate into a quark-antiquark pair, the second line gives the number of decay 
photons with an energy above 0.1 TeV, and the third line the number of photons above 0.2 TeV. 



Note the great sensitivity to the threshold : the number of photons above the treshold increases 
roughly as m 2 , and, as a result, the signal can only be detected if the WIMP mass is far above threshold. 
From Equation |?], and for our reference telescope, this signal is : 

Soft photon signal = / N 7 (E) dE = 8300 photons/yeara 2 / dn 7 (£) — T — ^ (8) 

JE threBh J 10 cm /s \ m x J 

The corresponding noise is the square root of the total number of background photons received above 
Ethresh ■ From Equations |l|, ^ and ||, it is approximately : 

f f m x ~) 0.2 TeV 

Soft photon noise = < / Nb ac k gr (E) dE > ~ 220 photons/ y'year (9) 

UE thre3h 'J Ethresh 

Note that we are conservative in integrating up to m x because the integral J dn 7 (E) is saturated well 
before m x , and therefore we are adding noise for no signal. 

The signal to noise ratio will be very low below m x ~ 1 TeV unless annihilation cross-sections are 
large. But they cannot be very far from < aV >haio— 10 -26 cm 3 /s , if WIMP's constitute the dark 
matter of the universe. WIMP's were produced in large numbers during the hot dense phase of the Big- 
Bang, and to constitute the dark matter, some fraction must have survived annihilation up to now. This 
annihilation becomes negligible when the temperature gets lower than a critical decoupling temperature, 
and the present relic density f2 x (in unit of the critical density) is j28| : 

_ , 2 2.5 x 10- 27 cm 3 /s , in . 
= — TTF^ (10) 

< (TV >BB 

up to logarithmic corrections (h is the Hubble constant divided by 100 km/s/Mpc, and we shall take 
h = 0.5). Note that the thermal average < aV >bb of the annihilation cross-section at the Big-Bang 
can differ from the thermal average < aV >haio in the halo, since the temperature at decoupling differs 
from the halo equivalent temperature. 

To account for dark matter in galactic halos, a mean density Qhaio — 0.1 is required [Q |j, ||], while 
a mean density £1 ~ 1 seems required on larger scales. When the WIMP annihilation rate is strong, 
il x << 1, and the density of WIMP's in the halo cannot be much larger than p x ~ Phaio^x/^haio ■ The 
remaining component could be due to brown dwarfs ||^] or to another WIMP. To interpolate the halo 
WIMP density between this strong annihilation rate and the weak annihilation rate (where p x ~ Phaio), 
Griest, Kamionkowski and Turner (GKT) (3(J suggested the formula p x ~ Phaio/[^ + ^haio/^x] > 
any smooth interpolation could be chosen. This GKT correction decreases the expected flux for a large 
annihilation cross-section, and Equation |^ becomes : 

Soft photon signal ~ 8300 photons/year a 2 f dn^(E) % ^ > haio f _V\ (11) 

J 'lO- 26 cm 3 /sV m x J [I + VLhaio/^x? 

"Line" signal 

Two WIMP's can directly annihilate into two photons (J dn 1 (E) = 2), and the energy of the photons 
is then equal to the mass m x of the (non-relativistic) WIMP, leading to a gamma line. A narrow line in 



5 



the range 10 2±1 GeV would provide a clear signature for dark matter annihilation (7[ ||, [| [l(], [H], |T^]. but 

the cross-section 027 for annihilation into 2 photons strongly depends on the WIMP model, and moreover 
it is expected to be much smaller than the annihilation cross-sections into a pair of quarks, leptons, Higgs 
or gauge bosons. For our reference telescope, the expected number of annihilation photons in the line at 

E = m x is : 

Line signals 166 photons/year a 2 > Ti° ( ) n i n 1 — 77712 ( 12 ) 

From Equation |[ the background noise for such a 7 line, with an energy bin AE = 3a = 0.3 E, is well 
approximated by : 

Line noise = U Nbackqr = 35 photons/^/year (13) 

m x 

Annihilations into a photon and another particle of mass M also gives a monochromatic photon, at an 
energy E ~ m x — M 2 /4m x . Some kind of spectroscopy is thus in principle possible, depending on the 



energy resolution. 



4 Neutralinos 

We now focus on supersymmetry E3 , the favorite extension of the standard model, because the lightest 
supersymmetric particle (LSP) X i s stable due to a conserved quantum number, R-parity, and because its 
relic density is naturally fl x ~ 1, making it a good dark matter candidate. There are two neutral gauge 
fermions B and W 3 and two neutral Higgs fermions H 1 and H 2 which mix when gauge symmetry is broken, 
and the four resulting mass eigenstates are called neutralinos. The lightest one, the neutralino, usually is 
the LSP. Its mass and couplings depend on 4 parameters : the masses Mi and M 2 of the gauge fermions 
B and W 3 before gauge symmetry breaking (usually related by the GUT relation Mi = §M 2 tan 2 9w)i 
the Higgs fermion mass [i and the ratio tan (3 of the vacuum expectation values of the two Higgs scalars 
Ju, 30, 3l], |32|, [H HJ]. In addition, the neutralino cross-sections depend on the masses of the quarks and 



leptons, and of their scalar partners, and on the Higgs bosons masses. There are weak direct experimental 
lower bounds from LEP on the neutralino mass fl35| but, assuming a GUT relation between the gluino 
mass M3 and the masses M\ and M 2 , more stringent indirect bounds can be obtained from UA2 |}6| and 
CDF |37[] experiments. It is then likely that the mass of the neutralino is larger than about 50 GeV. 

We are interested in the possibility that the neutralino be heavier than 100 to 200 GeV, the ACT 
threshold, because it allows a neutralino relic density near the critical density. We focus in this paper 
on the two extreme cases of a pure bino or a pure higgsino. A mixed state is not a priori excluded, but 
is quite unlikely in this high mass region : when ^> M 2 , the lightest neutralino is an almost pure B 
state, and an almost pure H when <?C Mi- The equations for annihilation cross-sections are extremely 
complex |}0[ [n], |32|, |33|, and we only give here a sketch of the results. 

Most of the time, a bino pair annihilates into a quark or a lepton pair, through the exchange of 
the corresponding s-quark or s-lepton of mass mj (m? > m x since, by assumption, \ is the LSP). The 
thermally averaged cross-section for binos at temperature T is approximately : 

<<rF>= c " 7 5 x 2V J 289-^ + 9160— ^ (14) 

216cos 4 0vf ( m f + m x> I m x m x \ m f + m x> J 

In this equation, a is the fine structure constant. As illustrated by Equation |lj, thermally averaged 
cross-sections for neutralino annihilation can be written as < aV >= a + bT/m x . Temperatures are very 
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low in the galactic halo and the a term dominates, whereas the b term dominates during the Big-Bang, 
down to the decoupling temperature T ( i ec — m x /2b. Equations [l^ and [l4| then yield a relation between 
the bino mass and its relic density Q x : 

m x = 0.56TcVh x /O x ~ (15) 



This value was obtained assuming that = m x , and it becomes much smaller if m~ t 3> m x . Therefore 
cosmology leads to an upper bound on the bino mass |30| [32], ||[ Q. It is very interesting that a 
similar upper bound is also indicated by naturalness arguments in the context of low energy supergravity 
models |3^] . A negligibly small relic density is excluded for a bino by the experimental lower bounds on 
the neutralino mass and binos should be around us in large numbers, even if they do not solve the dark 
matter problem. Assuming again the lightest possible s-top mass, namely = m x , the cross-section 
for annihilation in the halo into a top-antitop pair (for m x > m t0 p) is : 

<.V>_,2. 8 x 1 0-c ra »/,( T ^ Sv ) S (^X) 4 ( 16 ) 

if the bino is heavier than the top quark. The top quarks decay into a large number of soft photons, 
and protons, antiprotons, electrons, positrons, neutrinos and antineutrinos. The photon energy spectrum 
dn^/dE is computed by Lund Monte-Carlo J2lj. Figure 2 shows the expected excess in the number of 
photons in the direction of the center of the galaxy, seen by our reference telescope, as a function of the 
neutralino mass m x . The excess number of photons is compared to the background noise measured in 
directions away from the center of the Galaxy. The bino annihilation into soft photons will be difficult 
to detect, unless the dark matter density enhancement at the galactic center is larger than we assumed 
(i.e. if a > 1). 

Let us now turn to a H-like neutralino. The dominant higgsino annihilation channels are into WW 
or ZZ pairs, through the exchange of a chargino or a neutralino. The expression for its annihilation 
cross-sections is very complicated ||(], [H], |3^, |33j , but, in the mass range we are interested in, it scales 
roughly as : 

< aV >bb =s 2.2 x 1(T 26 cm 3 /s ( (17) 

leading to : 



rn v 



m x ~ STeVhy/Tl^ (18) 

The bound is larger than for the bino, because higgsinos are more strongly coupled than binos. Moreover, 
at low temperature, the higgsino cross-section scales as 1/m 2 (fermion exchange) while the bino cross- 
section scales as m 2 op /m 4 (scalar exchange). In the halo, the cross-section for higgsino pair annihilation 
into W's and Z's is : 

ttq 2 2 cos 4 6 W + 1 1 , 1Q , 
16 sin W cos 4 0w m x 

Higgsino annihilation is larger than bino annihilation, both in the Big-Bang (hence a smaller relic density 
for a given mass) and in the halo (hence a larger photon flux). Figure 2 indeed shows that higgsino 
annihilation into soft photons would be easier to detect than bino annihilation. However the signal is 
still marginal for our reference telescope. 

A gamma line would be a clearer signal . A neutralino pair cannot annihilate into two photons at the 
tree level, one loop graphs are required and probably lead to a tiny cross-section. There are many relevant 
graphs, in particular if annihilations into one photon and another particle (Z°, Higgs scalar) are included. 
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We are computing these graphs, but this long and tedious work is not yet completed. Nevertheless, an 
order of magnitude estimate may be derived from the graphs already computed p9[ . For the annihilation 
of a bino into two photons through quark/s-quark or lepton/s-lepton loops, this leads to : 

< a 2l V > halo ~ 4.2 x 10- 29 cm 3 /s ' (^) * ( 20 ) 

Of course, cancellations may happen between graphs, and the cross-section of Equation |2(] may be grossly 
overestimated. If it is not the case, Figure 3 shows that a bino could be detected at the 5 standard 
deviation level within one year of observation, down to the present 0.2 TeV threshold, again assuming 
TOf = m x . The cross-section for the annihilation of a higgsino into two photons is not yet known, and 
we boldly assumed the same cross-section as for a bino. Figure 3 shows that the GKT correction factor 
suppresses the higgsino line signal, due to the smaller higgsino relic density for a given mass. Thus a 
higgsino would be more difficult to detect in this way. The two searches, for a soft photon excess and 
for a line signal, are therefore complementary, as they are sensitive to different types of neutralinos, and 
they can be done simultaneously within the same experiment. 



5 Conclusions 

Dark matter made of particles in the 100 GeV - 1 TeV mass region could be detected through its 
annihilation products by an atmospheric Ccrcnkov telescope searching for energetic gamma rays. We 
apply this idea to the neutralino WIMP candidate from supersymmetry, and show that it may work, 
provided that : 

i) the neutralino is heavier than about 200 GeV, which is one of the windows allowing a relic density 
f2 x larger than 10~ 2 . 

ii) the density enhancement in the center of the Galaxy is as large as expected by the isothermal 
Ipser-Sikivie model. 

iii) a large enough atmospheric Cerenkov detector is built in the southern hemisphere, and aimed 
towards the galactic center. 

Improvements of the ACT technique to lower the energy threshold are of paramount importance to 
close the gap between the 0.05 TeV or so of present accelerator experiments, and the 0.2 TeV expected 
for our reference telescope. A complementary study could be based on the detection of the charged decay 
products, the electrons, protons and antiprotons. There is no directionality any more, and so no on-off 
experiment is possible, but consequently there is no prefered location for the apparatus. Integrating all 
over the galactic halo, preliminary calculations suggest that antiprotons (and protons) must display a 
strong enhancement below the neutralino mass, while other nuclei should not. 
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FIGURE CAPTIONS 

Figure 1 : Photon energy spectrum from the pair annihilation of neutralinos, showing the continuum 
due to quarks, leptons, Higgs and gauge boson decays, computed through a Lund Monte-Carlo, and the 
lines (convolved with experimental energy resolution) due to direct annihilations into two photons, and 
a heavy particle and a photon. The relative size of the lines and of the continuum is arbitrary. 

Figure 2 : Soft photon signal. Expected excess number of events per year above a 100 GeV threshold 
for a bino or a higgsino-like neutralino, as a function of the neutralino mass m x , compared to the noise 
(one-sigma fluctuations of the background). The bino and higgsino curves correspond to an Ipser-Sikivic 
model of the halo core with a = 1 , using the GKT correction to reduce the flux when the neutralino relic 
density Q x is smaller than the halo mean density Clhalo = 0.1 , and assuming a s-fermion mass m~ f = m x . 
The two vertical lines correspond to Q x h 2 — 1 for a bino and for a higgsino. 

Figure 3 : Line signal. Expected number of events per year in the line for a bino or higgsino-like 
neutralino, as a function of the neutralino mass m x , compared to the noise (one-sigma fluctuations of the 
background). The bino and higgsino curves correspond to an Ipser-Sikivie model of the halo core with 
a = 1, using the GKT correction, and assuming a s-fermion mass mj. = m x , and an energy resolution 
AE = 0.3 E. The two vertical lines correspond to VL x h 2 = 1 for a bino and for a higgsino. 
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